On the ground state and excitations of a four-component fermion model 



O 

o 

(N 



(N 
(N 



c3 



i 

c 

o 
o 



o 

o 
o 

03 



o 
o 



X 



You-Quan Li 1,2 , Guan-Shan Tian 1,3 , Michael Ma 1 ' 4 , and Hai-Qing Lin 1 

1 Department of physics, 
Chinese University of Hong Kong, Hong Kong, China 
2 Zhejiang Institute of Modern Physics, 
Zhejiang University, Hangzhou 310027, China 
3 Department of Physics, 
Peking University, Beijing, China 
4 Department of Physics, University of Cincinnati, 
Cincinnati, Ohio, 45221, USA 

(Dated: Received: February 2, 2008) 

The properties of the ground state and excitations of a four-component Hubbard-like model are 
studied by the Lieb-Schultz-Mattis approach. It is proven rigorously that the ground state of the 
model is nondegenerate and excitations are gapless at all band fillings. 
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The nature of the ground state and excitations of 
strongly correlated electronic systems plays an important 
role in our understanding of various fascinating phenom- 
ena such as metal-insulator transition, high-temperature 
superconductivity and colossal magnetoresistance. Re- 
cently, there has been much interests in the studies of 
correlated electrons with orbital degree of freedom^. In 
the case of double orbital degeneracy, a SU(4) theory was 
proposed as an ideal limit. [2j. 

Although there exists few rigorous results for general 
many-body Hamiltonians, significant progress has been 
made in one spatial dimension. In 1961, Lieb, Schultz, 
and Mattis (LSM) proved Q a remarkable theorem: For 
the spin 1/2 Hcisenberg antifcrromagnetic chain, the 
ground state is nondegenerate and the energy spectrum 
is gapless. This theorem has been extended recently to 
a hierarchy of generalized Heisenberg models 0, Q with 
high ranking symmetries, including the SU(4) model. It 
is of interest to obtain results for systems where the ki- 
netic terms caused by nearest neighbor hopping is not 
negligible. For the two-component Hubbard model, the 
non-degeneracy of the ground state was shown by Lieb|6J 
and gaplessness at 1/2 filling by Yamanaka, Oshikawa, 
and Affleck (YOA)0|. The extension of these results 
to four component Hubbard-type models is valuable be- 
cause such models not only describe electrons with dou- 
ble orbital degeneracy but also two layer systems with- 
out inter-layer hopping. Experimental realization of four 
component or two-band Hubbard model in one dimen- 
sion includes quasi one dimensional materials such as 
Na 2 Ti 2 Sb 2 amd NaV 2 O 5 0. 

In this work, we consider a four-component Hubbard 
model defined on a one-dimensional lattice with length L 

H = -t } j c+(x)c a (x') + U n a (x)n a >(x) 

. ° o<a' 

<x,x'> X 

+V n a (x)n a ,{x') , (1) 

<x,x'> 



where 



= 1,2,3,4; x identifies the lattice site, and 



(x,x') stands for nearest neighbor sites. c+(x) cre- 
ates a fermion of component a on site x, and n a {x) — 
c+(x)c a (x) is the corresponding number operator. The 
on-site interaction is assumed to be repulsive. The four 
components can be related to electrons with doubly or- 
bital degeneracy. Let us define the spin components by 
up (|) and down (J,) and the orbital components by top 
and bottom. Then the four possible electron singly occu- 



pied states arc |1) = | 
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Previous literature on such systems has largely con- 
centrated on the strong coupling limit of this Hamil- 
tonian, whereby it is mapped onto a 5(7(4) Heisen- 
berg mo del M, [lfll and related models away from SU(4) 
symmetry [ill Il2j. For the SU(4) Heisenberg Hamilto- 
nian, it has been established that there are gapless exci- 
tations at crystal momentum K = ±7r/2 and 7r[2j,|9j, 
In the case of no orbital degeneracy, it is well known that 
due to the existence of the charge gap for any U > 0, the 
low energy physics of the Hubbard model is essentially 
identical to the SU (2) Heisenberg model. It is of interest 
to investigate if this is the case for the SU (4) case also. 
In particular, much insight can be obtained if exact and 
rigorous results are available to complement approximate 
results using methods like bosonization |l3j . 

In the following, we first prove that the ground state of 
the Hamiltonian (JTJ with periodic boundary condition is 
nondegenerate when the total number of electrons is N — 
An with n being an odd integer. We then study whether 
the excitations are gapless for different band filling. 

To show nondegeneracy of the ground state, we 
shall apply a simplified version of Lieb and Mattis' 
method in proving the absence of ferromagnetism in one- 
dimensional itinerant electron lattice models [iH Hlj . 
First, we notice that, according to the representation 
theory of Lie algebra A%, for any multiplet, there is 
always one state lying in the subspace of zero weight 
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(0, 0, 0), which is spanned by the vectors with equal 
number of different kinds of fermions. Therefore, by in- 
troducing operators C|(x) = c\{x\)c\{x2) ■ ■ ■ c\{x n ) with 
x\ < X2 < ■ ■ ■ < x n , we can write the ground state wave- 
function as 



^C 1 T (x)Cl(x')C](x")Cl(x"')|0), (2) 



where /i stands for the set {x, x', x", x'"}. Obviously, ipo 
belongs to the subspace Vo(n), in which each vector has 
the same number n of a = 1, 2, 3 and 4 types of fermions. 
The total number of expansion terms in Eq. @ is (C£) 4 . 

In terms of this basis of vectors, we are able to write 
Hamiltonian into a matrix 7i which has the fol- 
lowing characteristics: (i) All of its off-diagonal ele- 
ments are either zero or — t, although its diagonal el- 
ements may have different signs. The reason why off- 
diagonal elements all have the same sign is due to the 
fact that, in Hamiltonian only hopping terms be- 
tween nearest-neighbor sites are present. In one dimen- 
sion, such hopping will not change the positional order 
of the fermions except at the boundary. However, if pe- 
riodic boundary condition is used for n odd, then the 
boundary hopping will not cause sign problems, (ii) Fur- 
thermore, TL is also irreducible. Namely, for any pair 
of indices m and n, there is a positive integer M such 
that the matrix element (7i M ) is nonzero since the 

V _ / run 

chain is connected by fermion hopping. Therefore, any 
pair of configurations C" 1 t (x)C 2 t (x')C ; |(x")C|(x / ")|0} and 
Cl(y)Cl(y')Cl(y")Cl(y"')\0) in subspace V (n) are con- 
nected by an appropriate number M of fermion hoppings. 

To such a Hermitian matrix, we are able to apply the 
well-known Perron- Frobenius theorem [lfij . It tells us 
that the expansion coefficients in Eq. 0) have the same 
sign and hence, ?/>o is nondegenerate in subspace Vo(n). 
To show that it is globally nondegenerate, we need to 
determine its quantum numbers by the continuity argu- 
ment. Adiabetically turning off the interactions of Hamil- 
tonian Q), we reduce it to a tight-binding Hq, whose 
ground state in Vo(n) is a spin-orbital singlet. It implies 
that tpo is also a spin-orbital singlet and hence, is globally 
nondegenerate . 

Next we study the nature of excitations for the case 
of N — 4n with n odd, focusing on whether they are 
gapless. Following LSM [g and YOA0, we consider the 
state 



\ip P ) = exp 



2tt \ ^ 
-l—r- y xPix) 



L 



x=\ 



\ip } = O(P)\ip ) 



(3) 



where operator P(x) = ^ Q v a n a (x), and the coefficients 
v a are to be chosen by the reasoning below. Our aim is 
to prove gaplessness. In order to do so, we want \il>p) 
to be orthogonal to \ipo) and to have an energy expecta- 
tion value which equals the ground state energy in the 
limit L — > oo. The operator O(P) "boosts" the crystal 
momentum of each a component particle by 2irv a /L. In 
order to preserve periodic boundary condition, v a has to 



be an integer. The state \ipp) then has a crystal momen- 
tum of K = ^ J2x=i P( x ) relative to that of the ground 
state. Since the Hamiltonian is translationally invariant, 
and the ground state is non-degenerate, it has a definite 
crystal momentum. Thus, a sufficient condition for the 
orthogonality condition is K ^ 2mir, where m is any in- 
teger, i.e., K should not be a reciprocal lattice vector. 
The expectation energy of \ipp) is given by 



where 



H' = HO 



Since the interacting part of H commutes with O the 
only difference between H' and H is in the tight binding 
part, and is given by a phase rotation of each c a {x) — > 
e~ l2lTVaX / L c a {x). This results in a change of the tight 
binding amplitude 



t -> t 



(«) 

x,x-\-l 



te 



-i27TV a I L 



x^L 



and t^\ = £ e l27 ™a( L - 1 )/ L _ Because v a is chosen to be an 
integer, periodic boundary condition is preserved, and 
this last bond will remain equivalent to all the other 
bonds. As a result, the difference H' — H will vanish 
as L — > oo. More precisely, using P{x) — ni(x) as an 
example, we have 



Ep-E a = -t(e 



-i2ir/L 



2n 



- l)^2(c\(x)ci(x + 1)) + c.c. 



= -2i(cos— - l)^2(c{(x)ci(x + 1)) 

x 

where the last line follows from {c\{x)c\(x + 1)} being 
real for a non-degenerate ground state. Actually, this 
method can only tell us that there is at least one low 
lying state with energy at most of order 1 / L higher than 
the ground state. Thus, the "gaplessness" may be due 
to having true gapless excitations or a discrete number 
of symmetry breaking states which are degenerate in the 
thermodynamic limit 0, . For this reason, our results 
below should all be taken implicitly as assuming there is 
no discrete symmetry breaking such as lattice translation. 

In particular we consider states generated by the fol- 
lowing four operators 



P^x) 

P2{X) 

Ps(x) 
Pi{x) 



= ni(x) , 

= m(x) + n 2 (x) , 

= ni(x) + n 2 (x) + n 3 (x) , 

= n 1 (x) + n 2 {x) + n 3 (x) + n 4 (x) 



(4) 



For general band filling, the crystal momentum K for 
all 4 states will not be equal to a reciprocal lattice vec- 
tor, and they are all orthogonal to the ground state, 
and there are gapless excitations at Ki generated by Pi. 
Band fillings equal to 1/4 (N — L), 1/2 (N = 2L), and 
3/4 (TV = 3L), however, must be discussed separately. 
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(i) 1/4-filled case: the ground state is a SU(4) sin- 
glet and we have ^ = N 2 = N 3 = N4 = N/4 = L/4, 
where N a = J2 X n a ( x ) ■ The crystal momenta of the 4 
states given by Eq. (0J are then tt/2, tt, Ztt/2, and 2tt 
respectively. Thus all 4 states are mutually orthogonal 
to each other, with the first three also orthogonal to the 
ground state. The fourth state cannot be concluded to 
be orthogonal to the ground state. Note that except for 
the non-interacting limit, the crystal momentum for each 
flavor is not a conserved quantity, only the total crystal 
momentum. Therefore, other combinations like 2n\ for 
example cannot be shown to be orthogonal to n\ + n-2- 

Our analysis establishes that at least at momenta tt/2, 
tt, 3tt/2 = —tt/2, there exists gapless excitations. In the 
large U limit, where the 4-component Hubbard model 
is mapped into an SU (4) symmetric Heisenberg model, 
it is known exactly that there are gapless excitations at 
K = 0, 7r/2, tt, and 3tt/2 and nowhere else. While our 
modified LMS approach cannot say anything about K = 
and cannot rule out gapless excitations at other K, it 
supports that the gapless modes are pinned at these K 
values for all U. 

(ii) 3/4-filled case: this case is equivalent to the 1/4- 
filled case by particle-hole symmetry. 

(iii) 1/2-filled case: the ground state is again a singlet 
and we have N t = N 2 = N 3 = N 4 = N/4 = L/2. The 
4 states from Eq. J3J now have momenta K = tt, 2tt, 37r, 
and 47T. Thus we can establish that gapless excitations 
must exist at K = tt. 

While gapless excitations are conclusively shown by 
the above analysis, they do not tell us much about the 
excitations beyond the crystal momentum. To try to 
understand these excitations more, we first consider the 
single band Hubbard model. By similar approach, one 
can show that twisted wavefunctions, for both spin up 
and down, 



one has 



IVO = ex P 



x = l 



(5) 



have crystal momentum K — tt at 1/2 filling and are 
orthogonal to the ground state. In the non-interacting 
limit, crystal momentum for each spin is conserved, and 
(ipn-f \ipni ) = O7 so there are two distinct gapless modes. 
In this limit, the gapless excitation with K = tt is a single 
particle-hole excitation of 2kp for one of the spins. These 
can then be cast in spin symmetric combination and spin 
antisymmetric combinations to give gapless charge and 
spin excitations. For non-zero U, the Hubbard model at 
1/2 filling always has a charge gap and only spin exci- 
tations are gapless. Correspondingly, only total crystal 
momentum is conserved and |VVi T ) and l^nt) no longer 
have to be orthogonal. Indeed, their overlap, (ipnt IV'ni ) 
is nonzero as verified by numerical calculations. Rewrite 

n T = -(n T + n L ) + -(n T - n L ) = \n + S z , 



exp 



2tt^ 



x=l 



exp 



exp 



-ij^xn(x) 



x=l 



x=l 



0{n)0{S z ) 



In the strong interaction limit, U — > 00, each site tends 
to be singly occupied, n{x) — > 1, so what left is the spin 
fluctuation S z (x). For finite U, it would seem that n\ 
generates both charge and spin excitations. However, in 
order for gaplessness, it must be that the overlap with 
charge excitations must vanish in the thermodynamics 
limit. While we cannot show this rigorously, we can get 
a partial understanding by considering the single mode 
approximation state 



|-01> 



L 

£ 

x = l 



exp (i7ra;)ni(a;)|^o) = 7;n(Tr)\ip ) + S z {tt)\%) 



So (^i) = I (n(-TT)n(Tr)) + (S z (-tt)S z (tt)) . Due to 



the charge gap, the density-density correlation is finite 
while the spin-spin correlation diverges as L — > 00. Thus 
there is no spectral weight for to be in a charge ex- 
citation. Conversely, our proof of gaplessness at tt shows 
that the spin-spin correlation function must diverge at tt. 



The twisted operator exp 



-i^£x=i*S»(*) 



is what 



Lieb, Schultz, and Mattis used in their paper (Ref. @) 
to show gapless excitation in the Heisenberg model. It 
cannot be shown to be orthogonal to the ground state in 
the Hubbard model because it has crystal momentum 0. 
Instead n a (x) must be used instead of S z (x). The fact 
that 0(n a ) gives a stronger result than that 0(n^ + njj 
and 0(ni — nj.) was noted by Yamanaka, Oshikawa, and 
Affleck 0. 

Returning to the 4-componcnt Hubbard model, we dis- 
cuss the analogous situation for the case of N = L/4. In 
this case, our proof shows gapless excitations for ±7r/2 
and tt. In the non-interacting limit, there are 4 orthogonal 
choices for P\ corresponding to whether we use n\nirij,, 
and 714, and they generate a single particle-hole excita- 
tion of 2kp. Linear combinations of these can be taken 
to form charge, spin, orbital, and spin-orbital excitations, 
the last three corresponding to the 3 diagonal generators 
of 517(4). For [7^0, the four n a states are not orthog- 
onal, and whether there are 3 or 4 gapless modes will 
depend on whether there is a charge gap (the other three 
are related by symmetry so they must be all gapless) . 
The "pinning" of gapless excitations at K = 2kp of the 
non-interacting system was used by YOA as a generalized 
definition of Luttinger's Theorem in ID. In analogy to the 
2-component Hubbard model, the gaplessness implies the 
spin-spin correlation function, the orbital-orbital correla- 
tion functions, and the spin-orbital-spin-orbital correla- 
tion functions must diverge at K = tt/2. In the strong 
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coupling limit, charge fluctuations are frozen, and we are 
definitely left with 3 gapless excitations, which are the 
three states with N x = N 2 = N 3 = N± = N/4 in the 15 
representation of SU(A). For K = it = 4kp, the gapless 
excitations for the non-interacting systems are pairs of 
2&f particle-hole excitations. The correlation functions 
that diverge at 4fcp are those corresponding to 4-particle 
Green's functions. 

For filling factor N a = p/q, with p and q integers, 
the orthogonality condition that K ^ 2mn will be sat- 
isfied for J2 a v a — 1>2, ...,q — 1. Accordingly there will 
be at least q momenta of gapless excitations. In the non- 
interacting limit, the J2 a v a — r gapless excitations are r 
particle-hole pairs. If we assume no broken translational 
invariance due to interactions (which should be the case 
for the Hubbard model with only on site repulsion) , then 
gapless excitations at these momenta persist when U 

In summary, we studied properties of the ground state 



and excitations of the four-component Hubbard model, 
in which electrons carry spin as well as orbital degrees of 
freedom. Considering cases where total number of elec- 
trons N — 4n with n odd for periodic boundary condi- 
tion, we showed rigorously that the ground state of this 
system is nondegenerate. Using the twist operators as 
specified in Eq. we addressed the issue of the exis- 
tence of gapless excitations. We showed that away from 
the filling factor 1/4, 1/2, and 3/4, the state produced by 
acting the twist operators (e.g., Eq. Q) on the nonde- 
generate ground state is orthogonal to the ground state 
and its variational energy approaches to the ground state 
energy in the thermodynamic limit. For the filling factor 
equals to 1/4 and 3/4, we showed that gapless excitations 
exist at crystal momenta tt/2, tt, 3tt/2 = — tt/2. 
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